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Abstract 

To estimate the deposition effect of PM2.5 (particle matter with aerodynamic diameter <2.5 |xm) in forests in northern 
China, we used the gradient method to measure the deposition velocity of PIV12.5 during the winter and spring above a 
deciduous forest in Olympic Forest Parl< and above a coniferous forest in Jiufeng National Forest Park. Six aerosol samplers 
were placed on two towers at each site at heights of 9, 12 and 15 m above the ground surface. The sample filters were 
exchanged every four hours at 6:00 AlVl, 10:00 AM, 2:00 PM, 6:00 PIVl, 10:00 PM, and 2:00 AlVl. The daytime and nighttime 
deposition velocities in Jiufeng Parl< and Olympic Park were compared in this study. The February deposition velocities in 
Jiufeng Park were 1.2±1.3 and 0.7±0.7 cm s~^ during the day and night, respectively. The IVlay deposition velocities in 
Olympic Park were 0.9±0.8 and 0.4±0.5 cm s~^ during the day and night, respectively. The May deposition velocities in 
Jiufeng Park were 1.1 ±1.2 and 0.6±0.5 cm s~^ during the day and night, respectively. The deposition velocities above 
Jiufeng National Forest Park were higher than those above Olympic Forest Park. The measured values were smaller than the 
simulated values obtained by the Ruijgrok et al. (1997) and Wesely et al. (1985) models. However, the reproducibility of the 
Ruijgrok et al. (1997) model was better than that of the Wesely et al. (1985) model. The Hicks et al. (1977) model was used 
to analyze additional forest parameters to calculate the PM2.5 deposition, which could better reflect the role of the forest in 
PM2.5 deposition. 
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Introduction 

Measuring the deposition effect of PM2.5 (particulate matter 
with aerodynamic diameter <2.5 |J,m) in forests is an important 
objective. This study is the first effort to estimate the deposition 
effect of PM2.5 in forests in northern China. Previous studies have 
been conducted in Europe and North America [1-2] but only one 
other study has been conducted in Asia. Kazuhide et al. [3] 
measured the deposition velocity of PM2.5 suffate in the summer 
using the gradient method in a deciduous forest at the eastern foot 
of Mt. Asama, Nagano Prefecture in central Japan. They obtained 
results that are similar to those measured in other parts of the 
world. Further investigations of deposition velocities in forests have 
been conducted in many sites in North America, Europe, 
Southeast Asia, and East Asia [4-6]. However, in East Asia 
deposition velocities in forests have only been studied in Japan and 
Chinese Taiwan locations [7-15]. 

The State Forestry Administration of China established this 
project to study the regulatory function and technology related to 
forest PM2.5. This study is a major project for exploring the role of 
the forest in PM2.5 deposition. 



Particulate matter 2.5 (PM2.5) is the most important contributor 
to haze. Under certain conditions, haze can cause the attenuation 
of atmospheric visibility. The use of motor vehicles in Beijing has 
significandy increased in recent years: since the 1 990s, the vehicle 
quantity has reached approximately 5.30 million [16]. Further- 
more, the total amount of atmospheric pollution in Beijing is 
increasing, which may be partiy due to the heavy use of 
firecrackers and fireworks on Chinese New Year and other similar 
occasions. Based on a report from the Ministry of Environmental 
Protection of China, the PM2.5 concentration reached 1000 ILg/ 
m' in Beijing on the Chinese New Year due to the use of 313,000 
boxes of fireworks [17]. The Chinese Government has been 
concerned about air pollution; thus, 35 environmental monitoring 
stations have been built in Beijing. These stations include 1 2 urban 
environmental monitoring stations, 1 1 suburban environmental 
monitoring stations, and 7 urban traffic environmental monitoring 
stations. A total of 24 hours of the measured data from all of these 
stations was published at http://zx.bjmemc.com.cn in 2013. In 
this study, we used the gradient method to survey temporal and 
spatial variations in PM2.5 concentration. 
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Experimental Section 

1. Ethics Statement 

This study did not involve any endangered or protected species. 
This work was conducted based on Forestry Standards "Obser- 
vation Methodology for Long-term Forest Ecosystem Research" of 
the People's Republic of China (LY/T 1952-2011). 

2. Sites 

The two field sites are depicted in Figure 1. One experiment 
site, Jiufeng National Forestry Park, is managed by the Forestry 
Committee of Beijing Forestry University and is available for 
teaching and research use by the university. It is located in Jiufeng 
National Forest Park in the Beijing Haidian District (40.06°N, 
116.09°E). It is a typical rural sampling site. 

The other experiment site, Olympic Forest Park, is managed by 
the Beijing Olympic Forest Park Development and Management 
Co., Ltd. It is located in a deciduous forest in the Olympic Forest 
Park in the Beijing Haidian District (40.258°N, 116.39°E). This 
site represents urban air pollution and is adjacent to Fifth Ring 
Road, one of the busiest roads in Beijing. Populus tomentosa is a 
winter-deciduous tree and is the dominant species in Olympic 
Forest Park. The dominant species in the Jiufeng National Forest 
Park is Platycladus orientalis evergreen. This site also contains Pinus 
tabulaeformis. The forests at both sites are classified as temperate 
forests, and flourish in similar climates. 

The height of the trees was approximately 8 m around the 
tower in both sites. The displacement height ranged approxi- 
mately from 4 m to 8 m (equations 1, 2 and 3). The roughness 
length ranged from 0.5 m to 1.5 m during the year for the 
Olympic and Jiufeng parks (equations 1, 2 and 3). According to 
measurements at the sites, the leaf area index [LAI] at the Jiufeng 
National Forest Park site was 3. 1 and 3.8 in the winter and spring, 
respectively. The leaf area index {LAI) at the Olympic Forest Park 
site was 1.8 in the spring. 



3. Sampling Program 

Two experiment stations were buHt, one in each forest, with 
iron towers that measured approximately 16 m in height. Figure 2 
depicts the setup of our experiment at both sites. 

PM2.5 fluxes were obtained using an atmospheric particulate 
sampler (KC-6120 integrated air sampler, Qingdao, Laoshan 
Electronic Instrument Factory Co., Ltd.). PM2.5 was collected on 
glass fiber filters (MK360, MunkteU&Filtrak GmbH, Sweden) at a 
flow rate of 100 L/min. The samplers were calibrated with the 
flow meter. Furthermore, the base of the filter film and the cutting 
head were ultrasonically cleaned with deionized water three times 
before each experiment. 

Three aerosol samplers were placed on each of the two towers at 
heights of 9, 12 and 15 m above the ground. The sample filters in 
the atmospheric particulate sampler were changed at 6:00 AM, 
10:00 AM, 2:00 PM, 6:00 PM, 10:00 PM, and 2:00 AM (i.e., 
three times during the day and three times during the night). The 
ultrasonic anemometer was placed on the iron tower at 1 5 m, and 
the meteorological instruments were placed on the iron tower at 9, 
1 2 and 15 m. The sampling times occurred in February 2013 and 
May 2013. The experiment was performed in the winter from 
February 22 to 28 and in the spring from May 7 to 12 in 2013. 
After every sample was obtained, each sample filter was sealed in a 
clean membrane polypropylene filter box to avoid contamination. 

4. IVleteorological Data and PM10 Data 

The meteorological data were measured by meteorological 
instruments in this study. The measured meteorological data 
included humidity and temperature (HMP45C, Campbell Scien- 
tific Inc., U.S.A.), wind speed and direction (014A/024A, Met 
One Instruments Inc., U.S.A.). An ultrasonic anemometer (Wind 
Master, GUI Instruments, Britain) was used to obtain friction 
velocity and the Monine-Obukhov length. 

Some data (PMIO and PM2.5 concentrations) were obtained 
from the Beijing municipal environmental monitoring center, the 
agency's official website is http:/ /zx.bjmemc.com.cn/. The data 




Figure 1. The map of the two fields. A: Olympic Forest Park in Beijing. B: Jiufeng National Forest Park in Beijing. 
doi:1 0.1 371/journal.pone.0097723.g001 
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Figure 2. The experiment setup at the study sites. The left panel is the tower in Olympic Forest Park in Beijing and the right panel is the tower 
in Jiufeng National Forest Park in Beijing. 
doi:1 0.1 371 /journal.pone.0097723.g002 



were applied to calculate the concentration ratios of PM2.5 from 
PMIO. Data is available in Data SI for this manuscript. 

5. Computation of PM2.5 Fluxes 

In this study, the gradient method was used to estimate PM2.5 
fluxes [18-20]. Table 1 lists all of the symbols and units used in 
this research. The equations 1, 2 and 3 were used to determine the 
d and Zo- 



used to determine the new value d, then the updated d value is 
substituted into the equations (2) and (3) to get a new g(d). We 
assumed the process has converged once the absolute value of the 
diflFerence between the updated g(d) and the previous one is less 
than 0.001. The final value oi d = g(d). 

The flux-gradient technique was used to determine the flux (F) 
from the measured vertical gradients of the concentration and the 
eddy diffusivity of sensible heat, as shown in equation (4). 



U{z) 



u" , Z — t 

— In 

k Zo 



(1) 



u(z) is the average wind velocity, u* is friction velocity, k is Von 
Karman constant, d is roughness length, ^5 is displacement height. 



fid\ - ^(^0-^(^2) _ HZi -d)-\n{Z2-d) 

■'^ ' C/(zi)-J7(z3) \n{Zi-d)-\n{Z3-d) ^ ' 

U(z\), U(Z2) and U(z:-,) are the wind velocity at heights Zi '^'^^ 
Zi, respectively. In this study we used 15, 12 and 9 m. 



(4) 



(5) 



g{d)=d- 



m 

fid) 



(3) 



f(d) and g(d) are arbitrary unction used in the iterative process and 
f(g) is the numerical derivative (the gradient between values in two 
iterations) of f(d). This equation set is used iteratively. First an 
initial value for d\s assumed. Then, the value oig(d) is updated and 



Ac = c(Z3)— c(Zi 



(6) 



In this equation, L is the Monin-Obukhov length, c is eddy 
concentration, Ac is the changes in the concentrations between ^1 
and Zi- is the integrated stability correction function in 
atmospheric deposition in relation to acidification and eutrophi- 
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Table 1. The parameters associated with the research. 



F 


fluxes 


Ac 


the changes in the concentrations between Zi and Z3 


u* 


friction velocity 


C 


the PM2.5 concentration at height Z3 


c* 


eddy concentration 


D 


the zero-plane displacement height 


LAI 


leaf area index 


L 


the Monine-Obukhov length 


k 


Von Karnnan constant 




the integrated stability correction function 


Z, 


height of 9 m 


D 


the transfer velocity 


Z2 


height of 12 m 


Zo 


the roughness length 


Z3 


height of 15 m 


RH 


relative humidity 


Vd 


deposition velocity 


Sc 


the Schmidt number 


P, 


the Prandtl number 


Rs 


the stomatal 


Rm 


the mesophyll 


Riu 


the outer surface resistances 


Rac 


resistances to transfer 




the resistance to transfer by buoyant convection 


Rgs 


resistances to uptake 


Rci 


the resistance to the uptake by exposed surfaces 


T 


temperature 







doi:l 0.1 371 /journal.pone.0097723.t001 
cation per Erisman and Draaijers (1995) [19]. 



From equations (4) and (5), F can be expressed by equation (9): 



= 21n 



l+.v 
2 



+ ln 



1-16 



1+^2 



[Z-d) 



- 2 arctan(x) + - (7) 



(8) 



-ku*Ac 



In 



(9) 



u* and L were averaged every 15 min, i*' was averaged every 4 h 
because Ac were from six periods: 6:00-10:00, 10:00-14:00, 
14:00-18:00, 18:00-22:00, 22:00-2:00, and 2:00-6:00. 

The deposition velocity, V^, was determined using the following 
equation adapted from Wesely and Hicks (1977) [21]: 
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Figure 3. The variations of temperature, relative humidity and wind velocity at the study sites. 

doi:1 0.1 371/journal.pone.0097723.g003 
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Vd-- 



F 
c 



(10) 



E 



z 



6. Empirical Models 

We chose three empirical models with which to compare our 
measured data: the models proposed by Wesely and Hicks (1977), 
Wesely et al. (1985), and Ruijgrok et al. (1997) [21,22,23]. The 
model by Wesely et al. (1985) is a dry deposition model over grass, 
whereas the Ruijgrok et al. (1997) model is a dry deposition model 
over forest canopy. The Wesely and Hicks (1977) model is a dry 
deposition model over all canopies. 

Slinn (1982) [24] and Wesely (1985) used the aerodynamic 
terms equation to express the deposition velocities of the aerosol 
particles. The following equation was used to calculate the 
deposition velocities, where is the aerodynamic drag, and V^^ 
is the surface deposition velocity: 



Vd-- 



1 



Rn 



(11) 



a. 
E 



E 



Furthermore, i?„ was calculated using the following equation, 
which is based on Erisman and Draaijers (1995) [25]. a* is the 
friction velocity, and If* is the integrated stability function for heat. 



+ 1 "N +1 



rM vo ^ 



^ m ro 



£ 

o 



In 



Rn 



ku* 



(12) 



Next, the Wesely and Hicks (1977) model was used for more 
complex parameters, where u is the wind speed, Sc is the Schmidt 
number, and Pr is the Prandd number: 



Vd-- 



1 



Ra + R,, + Rc 



(13) 



(5 
a. 



z 



Ru + Rh 



2 

ku* 



(14) 
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E 



1^ T— 



5 B 



2 a 

o 

<U _: 



1 

Rr 



1 



Rs+R.. 



1 

Rh, 



1 



1 



Rdc+Rcl Rac+R,s 



(15) 



The resistances in the upper canopy included R^ (the stomatal), 
R„ (the mesophyll), and Ri„ (the outer surface resistances). The 
resistances in the lower canopy were the resistance to transfer by 
buoyant convection {Rj^) and the resistance to the uptake by 
exposed surfaces [Rcij. The fourth term represents resistances to 
transfer (i?J and uptake (i?^„) on the ground [26]. 

The third value, Vj, was initially presented by Wesely et al. 
(1985) and was fit to the grassland ecosystem in the United States. 
This value can be calculated using equations 1 6 and 1 7 . The first 
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Figure 4. The mean PM2.5 concentration at tlie various KieigKits and sites. Tine mean PiV12.5 concentration at the heights of 9 m, 12 m, and 

1 5 m during the day and during the night in Jiufeng National Forest Park and Olympic Forest Park in February 201 3 and IVlay 201 3. The error bars are 
the standard errors at the various heights. (A: Olympic Forest Park in spring. B: Jiufeng National Forest Park in spring. C: Olympic Forest Park in winter 
D: Jiufeng National Forest Park in winter). 
doi:1 0.1 371 /journal.pone.0097723.g004 



equation was fit to the stability condition {L>0), and the second 
equation was fit to the instability condition {L<0). 



Vd-- 



w 
500 



(16) 



1+1 



300 



2/3 



(17) 




Jiufeng in spring Jiufeng in winter Olympic in spring Olympic in winter 



Different times and sites 

Figure 5. TKie concentration ratios of PM2.5 from PM10 at the study sites. 

doi:1 0.1 371/journal.pone.0097723.g005 
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Figure 6. The relationship between and V^. 

doi:1 0.1 371 /journal.pone.0097723.g006 



The third model was proposed by Ruijgrok and was fit to a 
European forest. If V, is the deposition velocity due to 
sedimentation, then V^^, is calculated using the following equation: 



Vd-- 



1 



+ Ra 



(18) 



RH>80,£ = 0.05m*" 



1 + 0.18EXP 



RH-SO 
20 



(21) 



RH< 80,^5 = 0.0067 



(22) 



Uh 



(19) 



RH>80,F, = 0.0067EXP('^^°^') (23) 



RH<80,E = 0.05i<*° 



(20) 



In the above equation, RH is the relative humidity, and u/, is the 
wind speed at the top of the canopy (approximately 9 m). 



Table 3. Mean deposition velocities, friction velocity and Monine-Obukhov length. 





Vd (cm s ^ 


u* (m s ^ 


L (m) 


Daytime in Olympic in Feb. 


0.02±0.03 


0.22±0.12 


-11.6±5.6 


Nighttime in Olympic in Feb. 


0.01 ±0.01 


0.12±0.06 


21.77±12.5 


Daytime in Jiufeng in Feb. 


1.2±1.3 


0.39±0.23 


-15.5±7.4 


Nighttime in Jiufeng in Feb. 


0.7±0.7 


0.29 ±0.1 3 


20.6±12.8 


Daytime in Olympic in May. 


0.9 ±0.8 ^^^^ 


0.35±0.22 


-17.15±9.7 


Nighttime in Olympic in May. 


0.4±0.5 


0.19±0.12 


27.97±14.3 


Daytime in Jiufeng in May. 


1.1 ±1.2 


0.42±0.25 


-15.9±7.8 


Nighttime in Jiufeng in May. 


0.6±0.5 


0.31 ±0.18 


22.8±6.4 


u* and L in the daytime and nighttime 


during the study period (February and May 2013). All of the 


values are the mean 


± standard deviation. 



doi:l 0.1 371/journal.pone.0097723.t003 
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Figure 7. The trend of the deposition velocities at the study sites and during the study periods compared with other models. A: 

Jiufeng Park in February. B: Jiufeng Park in iVlay. C: Olympic Park in May. 
doi:1 0.1 371 /journal.pone.0097723.g007 




Results and Discussion 

1. Atmospheric Conditions 

Meteorological conditions greatly impact the concentration of 
PM2.5 [27]. Therefore; therefore, the impacts of the conditions 
are important to analyze. Table 2 summarizes the meteorological 
conditions of Jiufeng National Forest Park and Olympic Forest 
Park, including temperature, humidity, wind speed, and solar 
radiation. Figure 3 presents the relationships between the 
meteorological conditions. We observed several common patterns 
between Jiufeng National Forest Park and Olympic Forest Park. 
The humidity and temperature were higher in the Olympic Forest 
Park, but the concentration of PM2.5, the solar radiation, and the 
wind speed were lower. 

2. PM2.5 Gradients 

Testing the changes in aerosol concentrations at multiple 
altitudes is difficult. The differences between the altitudes were 
minimal. A paired t-test reveals that the concentration at 1 5 m was 
significantly (P<0.01) higher than the concentration at 9 m during 
the spring (vegetation with young leaves), whereas the concentra- 
tions were not significantly different in the winter (vegetation 
without leaves). To decrease the sampling error, we used the same 
instrument to measure the concentration seven times at the same 
height of 9 m (or 15m). 

Fig. 4 shows the average concentrations measured at the various 
heights, sites and times. The concentration at 15 m was greater 
than that at 9 m, and concentrations appeared to decrease with 
height. This observation indicates a downward flux of the 
particulate aerosols (top down). The results indicate that the forest 
may have had a high absorption ability. As demonstrated in Fig. 4, 
during the daytime or nighttime, the average concentrations in 
Jiufeng National Forest Park in February significantly decreased 
with height, regardless of the time of day. This finding could be 



primarily attributable to the Platycladus orientalis, that is a type of 
evergreen species. In contrast, the gradients indicated no 
significant change during the daytime and nighttime in Olympic 
Forest Park, that could be primarily attributable to the presence 
and nature oi Populus tomentosa (a type of winter-deciduous tree). In 
conclusion, these results indicate that the evergreen forest had a 
higher absorption ability than the winter-deciduous forest in 
February. Similar results were obtained in May. These results may 
be because the Platycladus orientalis leaf structure is more complex 
than that of Populus tomentosa. The more complex the leaf structure, 
the more conducive the leaf is to adsorbing PM2.5. 

3. The Concentration Ratios of PIV12.5 from PM10 

The concentration ratios of PM2.5 and PM2.5-10 (PMIO 
minus PM2.5) to PMIO reflect different sources of pollution 
among regions and seasons. In the winter, natural gas heating is 
used throughout Beijing instead of coal; however, there is a lower 
socioeconomic belt around Beijing that uses old-fashioned 
fireplace burning coal for heating [28-29]. This population's 
reliance on coal heating is also responsible for the excessive sulfur 
dioxide discharge [30-31]. Another important source of pollution 
is the motor vehicle exhaust from the 5.3 million cars in Beijing 
[32]. 

As shown in Fig. 5, the concentration ratios of PM2.5 to PMIO 
were 42% in spring and 59% in winter in Jiufeng National Forest 
Park and 49% in the spring and 66% in the winter in Olympic 
Forest Park. This indicates that the environmental conditions had 
an enormous impact on the concentration ratios of PM2.5. 

4. The Change in the PM2.5 Concentration 

Fig. 4 illustrates the diurnal variation of the PM2.5 concentra- 
tions in various seasons. PM2.5 concentrations in the winter were 
much higher than those in the spring in Jiufeng National Forest 
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Park and Olympic Forest Park. The difference is due to the lack of 
leaves on the trees and the old-fashioned fireplaces (burning coal) 
in the areas circumjacent to Beijing, even though Beijing is a 
central heating city that uses gas [33]. Coal is the primary fuel for 
household heating and industrial production;burning coal pro- 
duces more sulfur dioxide than burning oil or gas. As shown in 
Figure 4, two diurnal peaks occurred in the PM2.5 concentrations 
at the test sites: 6:00-10:00 and 18:00-22:00. This phenomenon 
may be due to two major sources of PM2.5 those times. The major 
sources of PM2.5 were automobile exhaust and coal burning. 
Measurements indicated that during the morning and afternoon 
rush hours, the mean daily traffic flow of Fifth Ring Road in 
Beijing was 200 vehicles per minute; occasionally, the traffic jams 
occurred [34]. 

These values far exceed the PM2.5 standard set by the national 
AAQS of China [35]. PM2.5 emissions increased during heater 
usage and when meteorological conditions were unfavorable for 
atmospheric dispersion.In other studies [36-37], lower PM2.5 
concentrations were observed at special times, such as the 50th 
National Day in 1999, which can be attributed to specific 
procedures mandated by the government to reduce the emissions 
of PM2.5. Artificial simulation of rainfall was applied in the Che 
Gong Zhuang area on September 30, 1999. As a result, the PM2.5 
concentration was notably low (50 \lg/m'^). From October 3 to 12, 
the weather conditions were characterized either by still air or 
calm wind, and the PM2.5 concentration increased two-fold and 
reached approximately 160 Hg/m ' [36]. 

The PM2.5 concentrations gradually decreased from the winter 
to spring. The average PM2.5 concentration was 60% to 100% 
higher in th(" winter than in the spring, and 30% to 40% higher 
than the annual average [37] and was far greater than that of the 
EPA Ambient Air Quality Standard [38]. 

5. Deposition Velocities 

Table 3 lists the deposition velocities and relative parameters 
such as deposition velocity Vj (cm s~'), friction velocity u* (m s~'), 
and Monine-Obukhov length L (m), during the daytime and 
nighttime, in Olympic Forest Park and Jiufeng National Forest 
Park. A positive correlation exists between «* and V,j (Figure 6). 
The deposition velocities during the day and night in Olympic 
Forest Park were notably small in February because the vegetation 
in the forest lacked leaves. During February and May, the 
deposition velocities in Jiufeng were higher than in Olympic Forest 
Park. The main reason for this observation is that the PM2.5 
capture rate of needle-leaved evergreen forests is higher than that 
of broadleaved deciduous forests [39]. 

This study revealed that deposition increased in the daytime, 
whereas it decreased in the nighttime in Olympic Forest Park and 
Jiufeng National Forest Park. This result is consistent with previous 
studies. In Norway, Netherlands, America, Canada, Portugal [40], 
and Japan [10], the deposition velocities of aerosol particles were 
measured at the tops of coniferous forests. The measured values 
were greater than the computations by the Ruijgrok et al. (1997) 
model [23] in Kazuhide's study [3]. However, the deposition 
velocities observed in this study were comparable to those reported 
in Kazuhide's study [3]. Researchers from Japan suggested that 
th(' r<;as()n for this ambiguit)' could be a sampling error, as the 
measured values were larger than the modeled values. However, 
Horvath et al. [41] believed that it was important to correct the 
model parameters. Estimated Vj and deposition flux were strongly 
influenced by eddy diflFusivity in the roughness sub-layer [3]. The 
modification involves a height-dependent correction factor that 
ranges from 0.73 for ^= 22 m to 0.9 for ^— 34 m. In this study we 
used 0.64 [21,42] as a correction factor for the calculations. 



Therefore, the deposition flux increased 56% during the day and 
increased 52% during the night. The increase in deposition flux 
during the day was greater than that at night. The same general 
patterns were observed with the deposition velocity. 

6. Comparison between the Measured and 
Parameterized 

Because this type of data are not available for the Chinese 
mainland, it is very important to estimate the atmospheric 
deposition in this area. Figure 7. display the various deposition 
velocities at the study sites. The deposition velocities were 
calculated using the three models, and the error bars indicate 
the sampling errors in the experiment. We observed high 
deposition velocities during the daytime and low deposition 
velocities during the nighttime in Jiufeng National Forest Park 
and Olympic Forest Parks. The high deposition velocities in 
Jiufeng Park and the low deposition velocities in Olympic Park are 
consistent with the calculated results. However, the reproducibility 
of the Ruijgrok et al. (1997) model was relatively better than that 
of the Wesely et al. (1985) model. The calculated results from the 
Wesely et al. (1985) model were lower than those of the measured 
values. 

In the summer, Kazuhide et al. [3] measured the deposition 
velocity of the PM2.5 sulfate in a deciduous forest at the eastern 
foot of Mt. Asama, Nagano Prefecture, central Japan using the 
gradient method. Kazuhide et al. [3] explained the methodology 
and asserted that the botanical structure factor {LAI) was not an 
accurate parameterization. 

The measured value was higher than the calculated values from 
the Ruijgrok et al. (1997) model, the Hicks (1977) model and the 
Wesely et al. (1985) model. However, the Ruijgrok et al. (1997) 
model and the Hicks (1977) model were more suitable for Chinese 
forests and could better reflect the influence of the forest on PM2.5 
deposition. 

Conclusion 

The present study was conducted in Jiufeng National Forest 
Park and Olympic Forest Park in Beijing, China. To the best of 
our knowledge, this is the first report of the measurement of the 
deposition velocity of PM2.5 in China. This study also represents 
the first attempt to compare the deposition velocity in the city 
center with that in the suburbs. 

In general, the deposition was higher during the daytime than 
during the nighttime. likewise, the deposition was higher in the 
suburbs compared with the urban area. The deposition velocities 
of the aerosol particles were significantly higher in Jiufeng 
National Forest Park than in Olympic Forest Park during the 
same time periods. This may be due to the greater friction velocity 
in Jiufeng National Forest Park. Furthermore, the deposition 
velocities of the aerosol particles were also influenced by the eddy 
diffusion coefficient of the sub-layer surface roughness, as 
indicated by the report of Kazuhide et al [3] . 

Our results indicated that the deposition velocities of the aerosol 
particles were influenced by the friction velocity (Figure 6). The 
friction velocity was strongly influenced by the aerodynamic 
conditions [22,23,24,42]. Thus, the deposition velocities of the 
aerosol particles were strongly affected by aerodynamic conditions. 
These results are consistent with most recentiy reported results 
from other centers [22-24]. 

A potential sampling error resulted in arger measured values 
compared with the model values. Therefore, it was important to 
correct the model parameters to calculate the deposition velocities 
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of PM2.5 and to t:()mpare the calculated values of the models with 
the measured values. 

The measured results were more consistent with the Ruijgrok 

et al. (1997) model than with the Wesely et al. (1985) model. 
Different models were applicable to different profiles of the forests 
in different regions. The Ruijgrok et al. (1997) model was more 
applicable to the coniferous and broadleaved forests in northern 
China. 

Supporting Information 

Data SI Data SI is the data of meteorological data and 
PM2.5 concentration. This data contain seven parts: 1. 
meteorological datal include temperature, humidity, wind direc- 
tion,wind speed et al in different height in Olympic Forest Park. 2. 
meteorological data2 include temperature,humidity,wind direc- 
tion,wind speed et al in different height in Jiufeng National Forest 
Park. 3. met data include Monin-Obucov length and friction 
velocity. 4. PM2.5 DATAl is the PM2.5 concentradon in winter 
of Jiufeng National Forest Park. 5. PM2.0 DATA2 is the PM2.5 
concentration in spring of Olympic Forest Park. 6. PM2.5 DATA3 
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